Approved  for  public  release  by  the 
Director  of  the  National  Institute  of 
Standards  and  Technology  (NIST) 
on  October  9,  2015. 

NATIONAL  BUREAU  OF  STANDARDS  REPORT 

4795 


SECOND  EDITION  AUGUST  8,  1972 

SOME  OBSERVATIONS  ON  THE  USE  OF 
UNDERGROUND  RESERVOIRS  AS  HEAT  SINKS 


Report  to 


Office  of  the  Chief  of  Engineers 
Department  of  the  Army 
Washington,  D.  C.  20305 


U.S.  DEPARTMENT  OF  COMMERCE 

NATIONAL  BUREAU  OF  STANDARDS 


NATIONAL  BUREAU  OF  STANDARDS 


The  National  Bureau  of  Standards’  was  established  by  an  act  of  Congress  March  3, 
1901.  The  Bureau’s  overall  goal  is  to  strengthen  and  advance  the  Nation’s  science  and 
technology  and  facilitate  their  effective  application  for  public  benefit.  To  this  end,  the 
Bureau  conducts  research  and  provides:  (1)  a basis  for  the  Nation’s  physical  measure- 
ment system,  (2)  scientific  and  technological  services  for  industry  and  government,  (3) 
a technical  basis  for  equity  in  trade,  and  (4)  technical  services  to  promote  public  safety. 
The  Bureau  consists  of  the  Institute  for  Basic  Standards,  the  Institute  for  Materials 
Research,  the  Institute  for  Applied  Technology,  the  Center  for  Computer  Sciences  and 
Technology,  and  the  Office  for  Information  Programs. 

THE  INSTITUTE  FOR  BASIC  STANDARDS  provides  the  central  basis  within  the 
United  States  of  a complete  and  consistent  system  of  physical  measurement;  coordinates 
that  system  with  measurement  systems  of  other  nations;  and  furnishes  essential  services 
leading  to  accurate  and  uniform  physical  measurements  throughout  the  Nation’s  scien- 
tific community,  industry,  and  commerce.  The  Institute  consists  of  a Center  for  Radia- 
tion Research,  an  Office  of  Measurement  Services  and  the  following  divisions: 

Applied  Mathematics — Electricity — Heat — Mechanics — Optical  Physics — Linac 
Radiation1 2 — Nuclear  Radiation2 — Applied  Radiation2 — Quantum  Electronics3 — 
Electromagnetics3 — Time  and  Frequency3 — Laboratory  Astrophysics3 — Cryo- 
genics3. 

THE  INSTITUTE  FOR  MATERIALS  RESEARCH  conducts  materials  research  lead- 
ing to  improved  methods  of  measurement,  standards,  and  data  on  the  properties  of 
well-characterized  materials  needed  by  industry,  commerce,  educational  institutions,  and 
Government;  provides  advisory  and  research  services  to  other  Government  agencies; 
and  develops,  produces,  and  distributes  standard  reference  materials.  The  Institute  con- 
sists of  the  Office  of  Standard  Reference  Materials  and  the  following  divisions: 

Analytical  Chemistry — Polymers — Metallurgy — Inorganic  Materials — Reactor 
Radiation — Physical  Chemistry. 

THE  INSTITUTE  FOR  APPLIED  TECHNOLOGY  provides  technical  services  to  pro- 
mote the  use  of  available  technology  and  to  facilitate  technological  innovation  in  indus- 
try and  Government;  cooperates  with  public  and  private  organizations  leading  to  the 
development  of  technological  standards  (including  mandatory  safety  standards),  codes 
and  methods  of  test;  and  provides  technical  advice  and  services  to  Government  agencies 
upon  request.  The  Institute  also  monitors  NBS  engineering  standards  activities  and 
provides  liaison  between  NBS  and  national  and  international  engineering  standards 
bodies.  The  Institute  consists  of  the  following  technical  divisions  and  offices: 

Engineering  Standards  Services — Weights  and  Measures — Flammable  Fabrics — 
Invention  and  Innovation — Vehicle  Systems  Research — Product  Evaluation 
Technology — Building  Research — Electronic  Technology — Technical  Analysis — 
Measurement  Engineering. 

THE  CENTER  FOR  COMPUTER  SCIENCES  AND  TECHNOLOGY  conducts  re- 
search and  provides  technical  services  designed  to  aid  Government  agencies  in  improv- 
ing cost  effectiveness  in  the  conduct  of  their  programs  through  the  selection,  acquisition, 
and  effective  utilization  of  automatic  data  processing  equipment;  and  serves  as  the  prin- 
cipal focus  within  the  executive  branch  for  the  development  of  Federal  standards  for 
automatic  data  processing  equipment,  techniques,  and  computer  languages.  The  Center 
consists  of  the  following  offices  and  divisions: 

Information  Processing  Standards — Computer  Information — Computer  Services 
— Systems  Development — Information  Processing  Technology. 

THE  OFFICE  FOR  INFORMATION  PROGRAMS  promotes  optimum  dissemination 
and  accessibility  of  scientific  information  generated  within  NBS  and  other  agencies  of 
the  Federal  Government;  promotes  the  development  of  the  National  Standard  Reference 
Data  System  and  a system  of  information  analysis  centers  dealing  with  the  broader 
aspects  of  the  National  Measurement  System;  provides  appropriate  services  to  ensure 
that  the  NBS  staff  has  optimum  accessibility  to  the  scientific  information  of  the  world, 
and  directs  the  public  information  activities  of  the  Bureau.  The  Office  consists  of  the 
following  organizational  units: 

Office  of  Standard  Reference  Data — Office  of  Technical  Information  and 
Publications — Library — Office  of  Public  Information — Office  of  International 
Relations. 

1 Headquarters  and  Laboratories  at  Gaithersburg,  Maryland,  unless  otherwise  noted;  mailing  address  Washing- 
ton, D C.  20234. 

- Part  of  the  Center  for  Radiation  Research. 

3 Located  at  Boulder,  Colorado  80302. 


NOT  FOR  PUBLICATION  OR  FOR  REFERENCE 

NAHOiMAL,  BUREAU  OF  STANDARDS  REPORT 

4795 


SECOND  EDITION  AUGUST  8,  1S72 

SOME  OBSERVATIONS  ON  THE  USE  OF 
UNDERGROUND  RESERVOIRS  AS  HEAT  SINKS 


Report  to 


Office  of  the  Chief  of  Engineers 
Department  of  the  Army 
Washington,  D.  C.  20305 


U.S.  DEPARTMENT  OF  COMMERCE 

NATIONAL  BUREAU  Of  STANDARDS 


Preface 


This  report  is  being  reissued  at  this  time  because  of  the  current 
national  concern  for  energy  conservation.  The  technical  content  is 
pertinent  to  those  processes  where  available  thermal  energy  must  be 
stored  for  use  at  a later  time.  The  best  transfer  analysis  mathematic 
and  techniques  in  this  report  were  developed  in  regard  to  underground 
heat  sinks  but  can  be  used  for  problems  dealing  with  the  storage  of 
thermal  energy. 


. 

SOME  OBSERVATIONS  ON  TEE  USE  OP  UNDERGROUND 
RESERVOIRS  AS  ELAT  SINKS 

by 

B.  A,  Peavy  and  R.  S.  Dill 


Abstract 


Disposal  of  waste  heat  from  such  equipment  as 
Diesel  engines  and  refrigerating  or  air  condi- 
tioning machines  is  a practical  necessity  for 
the  occupation  of  large  underground  installa- 
tions and  the  operation  of  their  facilities. 
Cooling  water  supply  from  outside  may  be 
interrupted  during  emergencies  so  the  use  of 
underground  reservoirs  for  the  disposal  of 
waste  heat  Is  under  consideration.  Recircula- 
tion of  the  water  and  use  of  the  reservoir  as 
a heat  sink  have  the  advantages  that  the 
water  is  conserved  for  other  purposes  and  that 
the  total  heat  absorbing  capacity  is  augmented 
by  the  capacity  of  the  rock  surrounding  the 
reservoir  to  absorb  heat.  The  objectives  of 
the  present  work  were  (1)  to  find  whether  a 
reservoir  provided  for  an  existing  underground 
installation  is  large  enough  to  meet  design 
criteria  and  (2)  to  determine  for  the  general 
case  the  advantage,  in  terms  of  Increased  heat 
capacity,  of  recirculating  the  water  and 
utilizing  the  heat  absorbing  capacity  of  the 
surrounding  rock.  An  approximate  analytical 
solution  of  the  problem  was  sought  and  Its 
precision  was  investigated  by  means  of  tests 
with  the  large  reservoir  and  with  a smaller 
reservoir  provided  for  experimental  purpose.. 


INTRODUCTION 


This  paper  covers  analytical  and  experimental  work 
undertaken  to  establish  a means  of  predicting,  for 
erg  insuring  design  purposes,  the  heat  capacities  of 
v-ater -filled  underground  reservoirs.  The  capacity  of 
the  contained  water  is,  of  course,  the  product  of  its 
weight  and  temperature  rise,  but  the  total  heat  capacity 
of  the  reservoir  is  augmented  by  heat  transferred  to  the 
surrounding  rock  and  the  rate  of  heat  absorption  by  the 
surrounding  rock  is  dependent  upon  the  rate  of  tempera- 
ture rise  of  the  water  in  the  reservoir. 

It  is  probable  that  underground  reservoirs  will  be 
useful  as  emergency  heat  sinks  in  connection  with 
underground  protective  structures  assigned  to  resist 
attack  in  wartime,  since  heat  from  engines  or  from  air- 
conditioning  equipment  can  be  absorbed  by  such  a reser- 
voir if  the  outside  water  service  is  cut  off  due  to 
enemy  action  or  other  causes. 

An  analytical  method  for  predicting-  water  tempera- 
ture as  a function  of  time  and  heat  input  rate  was 
adapted  from  the  literature*  and  experiments  were  conducted 
with  an  existing  reservoir  to  corroborate  the  findings. 


MATHEMATICAL-  ANALYSIS  AND  ASSUMPTIONS 


To  develop  a theoretical  approach  to  this  problem, 
considerations  are  based  on  a hypothetical,  infinitely 
long  cylindrical  cavity  of  radius,  a,  from  which  heat  is 
transferred  radially  from  Its  lateral  surface  to  the 
surrounding  material  (rock)  that  is  assumed  to  be  bound- 
lees.  The  thermal  properties  of  the  surrounding  material 
o.re  thermal  conductivity,  k,  and  thermal  diffusi vi ty, ©o  . 

The  hollow  cylinder  contains  a fluid  (water)  of  specific 
heat,  c1,  Btu  per  pound  degree  F and  mass,  M1,  pounds 
per  foot  of  length  of  the  cylinder.  At  any  time,  t,  the 
temperature  6 of  the  fluid  is  assumed  to  be  the  same  as 
the  temperature  at  the  surface  of  the  cylinder.  A con- 
stant heat  flux  Q.,  Btu  per  hour  per  foot  length  of  cylinder, 
is  added  to  the  fluid  from  ah  external  source  and  the  heat 
balance  becomes; 
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From  Carslaw  and  Jaeger,  “Conduction  of  Heat  in  Solids", 
p.  2 dii , the  temperature  of  the  fluid  at  any  time,  t,  is; 
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Symbols  appear  at  the  end  of  this  report. 

Tunnels  used  for  underground  reservoirs  have  irregular 
or  non-circular  cross  sections,  and  finite  lengths . The 
equivalent  radius  of  the  tunnel  is  computed  from  che 
average  perimeter  of  the  actual  irregular  cross  section, 
hxcept  for  cases  in  wnich  the  width  to  height  ratio  is 
large,  this  gives  approximately  the  same  heat  flow  path 
as  the  cylindrical  case.  For  a reservoir  or  tunnel  of 
finite  length  some  heat  is  transferred  into  the  rock  at 
the  closed  end  and  into  the  rock  surrounding  the  opening 
at  the  open  end.  This  heat  is  not  accounted  for  in  the 
theoretical  treatment  and  can  often  be  neglecced  because 
(1)  this  heat  provides  a factor  of  safety  in  design 
computations,  (2)  this  heat  is  negligible  \%rhen  the  ratios 
of  length  to  width  or  height  are  great,  since  in  such  a 


ca-o  the  surface  area  of  tuo  cna  is  small  in  comparison 
zo  the  lateral  surface  area,  and  (3)  relative  time  for 
use  of  the  reservoir  is  small. 


EXPiR I MENTAL  UNDERGROUND-  hEBERVOIR 


The  experimental  underground  reservoir  is  66  foot 
long  with  an  irregular  cross  section  on  the  order  of 
seven  feet  square.  The  equivalent  cylindrical  radius 
was  computed  to  be  I4.2  feet  from  averaging  perimeter 
measurements  taken  every  two  feet- along  the  length. 

The  reservoir  is  surrounded  by  greenstone  except  for 
one  end  which  was  closed  during  tests  by  a three-foot 
concrete  dam  with  a plywood  diaphragm  above  it.  The 
properties  of  the  greenstone  are:  thermal  conductivity, 

1.1-5  Btu/hr-f  t-{  deg  F);  specific  heat,  0.2  Btu/lb 
(aeg  F);  ana  density,  I06  lbs/ft3.  Water  at  the  dam 
was  pumped  through  a tank  containing  electric  heaters 
and  discharged  back  to  the  reservoir.  Pleat  input  to 
the  water  by  electric 'heater s and  pump  was  measured  by 
a watthour  meter.  Temperature  measurements  were  made 
at  various  positions  in  the  rock,  at  the  rock  surface 
and  in  the  water  by  means  of  thermocouples. 


OBSERVATIONS  WITH  FOG-NOZZLE  GRAYING  OF  HOCK 
SURFACES  ABOVE  WATER  LEVEL 


For  the  duration  of  this  test,  approximately 
11,390  gallons  or  95*000  pounds  of  water  was  maintained 
in  the  reservoir.  The  water  was  discharged  back  to  the 
reservoir  through  a two-inch  header  with  fog-nozzles 
spaced  every  seven  feet  along  the  length  of  the  reser- 
voir. Nozzles  were  directed  at  the  exposed  rock 
surface  above  the  water  level. 

Heat  was  added  to  the  water  at  a rate  of  55*^00 
Btu  per  hour  for  hours.  Figure  2 shows  the  ime- 

temperature  relationship  for  the  water,  rock  surface 
and  rock  at  depths  of  2 and  5 feet.  The  values  for 
the  dimensionless  groups  of  equation  (1)  are: 
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The  result  of  substitution  of  experimental  data  in  the 
above  equations  is  shown  on  figure  3 along  with  the 
theoretical  plot  for  G=3«5* 

The  experimental  results  agree  closely  with  those 
yielded  by  equation  (1)  closely  up  to  ect/a^  = 0.2. 

The  discrepancy  between  the  theoretical  and  experimental 
results  for  values  of  ^ct/a2  greater  than  0.2,  can  be 
explained  by  the  fact  that  the  theoretical  treatment 
does  not  take  into  account  the  heat  transferred  into  the 
end  of  the  reservoir.  This  can  be  demonstrated  more 
clearly  by  comparing  the  heat  balances  as  follows: 

Heat  Balance  Based  on  At  End  of  250  Hours 

Temperature  Increases 
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Theoretical  Heat  Balance  Based  on  liquation  (1)  At 
End  ox  2$S  Hours 

heat  in  water  3*45  * 10^  Btu 

L 

heat  in  rock  cylinder  10.75  * 10  Btu 

Total  ll|.2  x 106  Btu 

The  heat  balances  show  that  for  the  same  amount  of 
heat  added  to  the  system,  the  distribution  assumed  for 
the  theoretical  case  was  different  from  that  which 
actually  occurred  during  the  experiment  and  that  the 
water  temperature  rise  was  less  for  the  experimental 
case.  On  the  basis  of  the  total  amount  of  heat  added, 
the  theoretical  heat  absorbed  by  the  water  was  approxi- 
mately 3 percent  greater  than  the  experimental.  An 
empirically  determined  relationship  for  values  of 
*ct/a^  greater  than  0.2  can  be  written: 

/0k}  ™ / Qk  n 0 s Q3kt  0 > 

^ 'actual  " 'theoretical  ” M * ^ 


which  closely  follows  the  experimental  plot  beyond 
«s:  t / a2  — 0.2. 

Unaccountable  heat  losses  from  the  system  include 
water  seepage  into  the  reservoir  from  fissures  in  the 
rock  and  out  the  overflow  pipe,  and  water  vapor  transfer 
and  heat  loss  through  openings  above  the  dam.  As  a 
percentage  of  the  total  heat  added  these  losses  can  be 
considered  negligible. 

Q3SERVATI0 N S WITHOUT  3 PRAYING  OF  SURFACES 
ABOVE  WATER  LEVEL 


For  the  duration  of  this  test,  approximately 
15,200  gallons  or  120,600  pounds  of  water  was  maintain©- 
in  the  reservoir-.  Circulating  water  was  drawn  from  the 
reservoir  near'  the  bottom  of  the  dam,  passed  through 
the  heater  and  pump  and  returned  to  the  reservoir 
through  a long  pipe  which  extended  the  length  of  the 
tunnel  and  delivered  the  water  at  the  far  ena. 
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Heat  was  added  to  the  circulating  water  at  a T ate 
of  61; , 500  3tu/nr  for  256  hours  and  figure  6 shows  the 
time- temperature  relationship  obtained  for  the  water, 
rock  surfaces  arid  rock  depths  of  2 and  5 feet.  The 
values  for  the  dimensionless  groups  of  equation  (1) 
ore: 
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The  result  of  substitution  of  experimental  data  in  the 
above  equations  is  shown  on  figure  5 along  with  the 
theoretical  plots  for  G=2 . 5 and  3.0.  The  experimental 
results  are  slightly  above  those  yielded  by  equation 
(1)  for  G=2,76  and  for  values  of  «ct/a^  less  than  0.6. 

For  greater  values  of  ect/a^  they  fail  below  the  results 
of  equation  (1).  The  discrepancy  between  the  theoretical 
and  experimental  results  when  <ct/a^  exceeds  0.6  can  be 
reconciled  as  in  the  case  of  the  spraying  test  by  the 
following  empirical  relationship: 


/■  ©k)  = / ©k y __  0. 02kt 

q ;actual  Q,  theoretical  M"7 


DISCUSSION  AND  CONCLUSIONS 


Two  assumptions  are  essential  for  this  analytical 
study  of  the  problem  which  cannot  be  met  in  a practical 
sense.  They  are  as  follows: 
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1.  That  heat  is  transf erred  radially  from  the 
lateral  surface  of  the  hollow  cylinder  and  that  the 
ends  of  the  hollow  cylinder  are  insulated  from  heat 
flow.  This  was  not  true  during  the  experiments 
because  there  was  heat  flow  into  the  ends  and  its 
effect  gave  lower  water  temperatures  than  predicted 
by  theory. 

2.  That  water  is  so  well  stirred  that  the  tempera- 
ture of  the  water  is  the  same  as  the  temperature  of  the 
surface  of  the  reservoir.  This  impl5.es  that  there  is 

no  thermal  resistance  between  the  water  and  the  rock 
surface.  In  the  two  tests  (figures  2 and  I4 } there  was 
a temperature  difference  between  the  water  and  the 
surface.  This  effect  vias  more  evident  when  there  was 
no  spraying  of  rock  surfaces  above  the  water  level. 

The  above  two  effects  appear  to  counteract  each 
other  and  the  greatest  effect  after  a period  of  time 
was  shown  to  .be  heat  stored  in  the  end  of  the  reservoir. 
This  was  borne  out  by  lower  water  temperatures  than 
predicted  by  the  theoretical  approach,  as  shown  in  the 
empirical  equations  (2)  and  (3). 


COMPARISON  OF  RESULTS  WITH  AND  WITHOUT  SPRAYING 

OF  ROCK  SURFACES 


Comparison  between  the  two  tests  cannot  be  directly 
made  because  of  the  difference  in  heat  input  rates,  but 
when  equation  (2)  is  subtracted  from  (3),  the  following 
relation  results: 


0 , OlQt 
M 1 


where  9.  and  9^  are  the  temperature  rises  fer  the  non- 
spraying  and  spraying  tests,  respectively.  By  assuming 
a steady  heat  flux,  Q,  of  700  Btu/hr-ft,  mass,  M 1 , of 
1350  lbs /ft,  and  time,  t,  of  ipQG  hour  3 5 the  temperature 
difference  between  the  two  cases  becomes  2F.  This 
represents  an  increase  of  approximately  16  hours 
attributable  to  spraying  water  on  rock  surfaces  above  the 
water  level. 


DLitJIGN  APPLICATIONS 


The  results  of  the  experiments  indicate  that 
equation  (1)  and  figure  1 can  be  used  for  estimating 
heat  capacities  of  underground  reservoirs  with  reason- 
able accuracy  in  most  practical  cases.  For  values  of 
exceeding  0.2,  the  Pieat  capacity  of  the  rock 
surrounding  the  ends  of  a reservoir  becomes  appreciable 
and  makes  the  total  capacity  greater  than  that  predicted 
by  equation  (1)  and  figure  1.  Correction  factors  are 
suggested  in  equations (2)  ana  (3).  However,  under- 
ground reservoirs  used  as  heat  sinks  are  likely  to  be 
larger  -ihsun  the  experimental  reservoirs  so  that  values 
of  exceeding  0.2  are  not  likely  in  practice  and 

the  correction  is  expected  to  be  unnecessary  in  most 
cases.  For  the  sake  of  example,  assume  an  equivalent 
radius  of  10  feet  and  «?=  O.Ob  ft^/br,  men  the  time 
necessary  to  reach  a value  of  &?t/a^  = 0.2  is  500  hours 
or  21  days.  This  is  probably  a considerably  longer 
time  than  is  anticipated  for  the  use  of  a reservoir 
a.a  an  emergency  heat  sink. 

F.or  a sample  problem,  consider  an  underground 
reservoir  with  dimensions  230'  long  by  2k  1 wide  and 
29*  high  holding  1.5  million  gallons  of  water.  The 
thermal  properties  of  the  rock  are  the  same  as  that 
for  the  experiments  and  the  initial  water  temperature 
is  5 OF.  The  value  for  the  equivalent  radius  is; 
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input  rate  to  the  water  of  6 x 10°  Btu/hr 
the  dimensionless  groups  are: 


th 


2 Tf  (20)2  x l0l\5  x 2 30 

1,500,000  X 0.33  X 1.0  X 0.0395 


x 


7 


4 *-3T>  to 


-10 


oC  t 


0 . 0 3 9 9 1 

~T2oT^~ 


= x io*^t 


Gk  90  x 1,45.x  230 

§ = 6 i "lO5 


0.00278 


kef erring  to  figure  1A,  cct/a^  = 0.012,  and  the  tirae 
equals  121  hours  or  a savings  of  16  hours  due  to  heat  transfer 
to  the  surrounding  rock. 

The  above  example  is  taken  from  an  actual  under- 
ground reservoir  where  results  of  a test  on  this 
reservoir  closely  substantiate  the  above  figures. 

The  value  of  G with 'the  reservoir  completely 
filled  is: 


G = S££- 

u p^c 5 


with  water  in  the  reservoir 


pc 

G ~ 31*2 

where  p and  c are  the  density  and  specific  heat  of  the 
urrounding  rock,  respectively.  For  these  experiments 
= 186  lbs/ft-3  and  c = 0.2  Stu/lb-°F,  then  G = 1.2. 
her ef ore,  optimum  conditions  are  reached  when  the 
eservoir  is  completely  filled  with  water. 
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SYMBOL  6 

Symbols  and  constants  used  for  the  underground 
reservoir  are  as  follows: 

a = radius  of  equivalent  cylinder,  feet 
= I4.2  feet 

c»  = specific  heat  of  fluid,  Btu/lb  - (deg  F) 

= 1.0  (water) 

G = mathematical  quantity,  used  in  equation  (1) 

k = thermal  conductivity  of  rock,  Btu/nr  - ft 

- (deg  F) 

= 1.1*5  (greenstone) 

M*  = weight  of  fluid  per  unit  length  of  cylinder 

WnTt. 

Q,  = heat  flux  to  fluid  from  external  source, 
Btu/hr  - ft 

t = time,  hrs 

p 

**  = thermal  diffusivity  of  rock,  ft  “/hr 

= 0.0395  for  greenstone 

pc  = volume  specific  heat  of  rock,  Btu/ft^  - 
(deg  ?) 

= 36.7  (greenstone) 

p!c*  = volume  specific  heat  of  fluid,  Btu/f tJ  - 
(deg  F) 

= 62. ^ (water) 


temperature  rise  of  water  above  the  initial 
temperature  at  time,  t,  deg  F. 
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